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I— I Abstract 

p ^ In this paper we provide evidence that financial option markets for equity indices 

give rise to non-trivial dependency structures between its constituents. Thus, if the in- 
dividual constituent distributions of an equity index are inferred from the single-stock 
option markets and combined via a Gaussian copula, for example, one fails to explain 
the steepness of the observed volatility skew of the index. Intuitively, index option 
prices are encoding higher correlations in cases where the option is particularly sensi- 
j>. tive to stress scenarios of the market. As a result, more complex dependency structures 

^-H emerge than the ones described by Gaussian copulas or (state-independent) linear cor- 

relation structures. In this paper we "decode" the index option market and extract 
this correlation information in order to extend the multi-asset version of Dupire's "lo- 
cal volatility" model by making correlations a dynamic variable of the market. A 
^> "local correlation" model (LCM) is introduced for the pricing of multi-asset deriva- 

fives. LCM achieves consistency with both the constituent- and index option markets 
• . by construction while preserving the efficiency and easy implementation of Dupire's 

. ^ model. 

^ Keywords: implied correlation, local correlation, stochastic correlation, correlation 

skew, index skew, basket options, multi-asset Dupire model, multi-asset local vol 
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1 Introduction 



In 2008, the crisis in credit markets led to a spike in correlations that resulted in large losses 
for several single-stock exotic trading books. As most market participants did not anticipate 
the crises, one is tempted to also put these equity losses under the category of unpredictable 
and unforeseen "bad" events. We do not fully subscribe to this view. 

In this paper we argue that, at least to some extent, the losses could have been allevi- 
ated if models had incorporated non-trivial correlation information, already encoded in the 
index option market, properly into the pricing of multi-asset derivatives. This is because the 
index volatility market already implies rising correlations in downturn scenarios. 

Most investment banks utilise a Gaussian copula for the pricing of multi-asset exotic 
European payoffs. For more involved payoffs a multi-asset version of Dupire's local volatil- 
ity model has become the market standard for the pricing and hedging of those products. 
The model infers a local spot-dependent volatility surface from option prices for each as- 
set and combines them via a (constant) correlation matrix between its Wiener increments. 
However even when correlations are assumed time-dependent, the model does not incorpo- 
rate all dependency information inherent in the derivatives market properly. In particular 
with downward scenarios, index options imply correlations that can be substantially higher 
compared to the "at-the-money" case. As a result the cost of short cross-gamma position 
hedging could be priced in advance, more appropriately, by means of a better model. 

In chapter 2 of this paper we take a snapshot of the options market and deduce the 
distribution of the Dow Jones Euro Stoxx 50 (.STOXX50E) and DAX (.GDAXI) indices 
together with all the composite distributions. In both cases we show that only half of the 
volatility skew of the index is explained by the skew of the index constituents. We conclude 
that (term) correlation effectively depends on the strike of the option which substantially 
impacts the distribution of the index on the downside. 

It is important that multi-asset derivative models incorporate these effects even if a 
product is only sensitive to a subset of stocks in the index. This is because the picture 
remains fundamentally the same: stocks correlate more strongly in a bear market, when 
investment behaviour is driven by fear, than during bullish times. As a result portfolio 
managers may find themselves less diversified, even when the composition of their portfolio 
stays the same. As far as hedging of derivatives is concerned, an increase in correlations 
has the unpleasant effect of increasing the cost of managing short cross-gamma positions. 
This arises as a result of the non-linear behaviour of (multi-asset) options' prices and man- 
ifests itself in the fact that the option's sensitivity to one asset increases with the value 
of a second decreasing. During the credit crunch for example, the Nikkei often finished 
down for the day creating a long delta position in EuroStoxx 50 for a hedging book in 
London, even though the market there was closed. As a result, the trader would need 
to immediately sell EuroStoxx 50 at the open. However EuroStoxx 50 would be hkely to 
open down as well because of the large positive correlation between the two assets. This 
gives rise to systematic losses during periods when the correlation exceeds expectations. 

The goal of this paper is to develop a model that incorporates higher correlations in down- 
ward scenarios by "decoding" the state-dependent correlation information of the index option 
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market. One expects that such a model anticipates losses such as those described above at 
the inception of the trade and adjusts the pricing accordingly. In addition one would expect 
more accurate deltas in this case. 

In order to make LCM a model that can be useful for practical applications, one needs 
to provide a robust as well as efficient implementation for the underlying algorithms. This 
is the strenght of Dupire's local volatility model which is simple to implement, numerically 
sufficiently efficient and provides consistency with the volatility skew of individual option 
markets at the same time. The goal of this paper is to extend this methodology to the 
multi-asset case by providing consistency with the index option market in addition to the 
individual option markets. Note that, similar to the Dupire framework, LCM provides a 
particular way of achieving consistency with the option markets, however this solution is 
by no means unique. Chapter 7 provides alternative solutions to the same problem which 
allows one to quantify "residual exotic correlation risk"'. 

Work in this direction has also been pursued by others. See for example [1] , [2] , [2] , [1] , [5] . 

Chapter 3 lays down the ground work for the "local correlation" model (LCM) that is 
introduced in Chapter 4. A parameterisation of correlations is introduced that preserves the 
positive semi-definiteness and displays some notion of strictly increasing correlation matrices. 
Chapter 5 provides simulation results of LCM and provides explicit results for the strike 
dependence of correlation. Chapter 6 deals with some implementation details and chapter 7 
discusses residual correlation risk. Chapter 8 concudes our analysis. 
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Ticker Wgt 


Ticker Wgt 


Ticker Wgt 


Ticker Wgt 


Ticker Wgt 


AEGN.AS 2.67 

GASI.MI 2.26 
BAYG.DE 1.57 
DAIGn.DE 1.69 

ENEI.MI 12.11 
DANO.PA 1.14 
OREP.PA 0.45 

REP.MC 1.50 
SIEGn.DE 1.55 

CRDI.MI 30.25 


AIRP.PA 0.50 
AXAF.PA 3.17 
BNPP.PA 1.66 
DBKGn.DE 1.18 
ENI.MI 4.63 
IBE.MC 6.54 
LVMH.PA 0.49 
RWEG.DE 0.77 
SOGN.PA 1.10 
UNc.AS 2.95 


ALVG.DE 0.86 
BBVA.MC 7.12 
CARR.PA 1.16 
DBlGn.DE 0.37 
FOR.BR 4.16 
ING.AS 3.95 
MUVGn.DE 0.38 
SASY.PA 1.99 
TLIT.MI 19.43 
SGEF.PA 0.94 


ALSO. PA 0.38 
SAN.MC 15.50 
SGOB.PA 0.76 
DTEGn.DE 5.66 
FTE.PA 3.61 
ISP.MI 18.33 
NOKIV.HE 7.12 
SAPG.DE 1.66 
TEF.MC 7.79 
VIV.PA 2.33 


ISPA.AS 1.76 
BASF.DE 1.75 
CAGR.PA 1.93 
EONGn.DE 3.80 
GSZ.PA 2.53 
PHG.AS 1.85 
RENA.PA 0.38 
SCHN.PA 0.49 
TOTF.PA 3.99 
VOWG.DE 0.16 



Table 1: Decomposition of EuroStoxx 50 on the 31 July 2009. The first column displays the 
(Reuters) names of the assets followed by the weight in the index. 



2 "Decoding" the Index Options Market 

In this chapter we take a snapshot of the market as of 31 July 2009 and try to explain 
the observed index option prices of EuroStoxx 50 and DAX in terms of their constituent 
distributions. One finds that the skewness of the individual distributions fails to fully explain 
the volatility skew of the index basket itself. Hence, the index volatility market encodes non- 
trivial dependencies between its constituents. Let It denote the value of the index at time 
T. The price of a call option on the index with maturity T and strike K is given by C(T, K). 
Note that 



/>oo 

C{T,K) = DF E^{{It- K)+) = DF dl prob{T, - (1) 

Jo 

where DF is the discount factor to maturity and denotes the expectation value under 
the forward adjusted measure Q. Thus, Eq|l] relates option prices to the (forward-adjusted) 
implied density prob{T, It) that is inferred by the option market with maturity T. Assume 
that for a given maturity T all strikes are traded. In this case Eqjl] can easily be inverted 
with the result 

' B'^^T'^K,roHT,K). (2) 



DF dK^ 



Let us denote by Si{T), i = 1, .., n, the price of constituent i at time T. The index is generally 
defined as a weighted sum of its constituents e.g. 



It = J2^^S,{t). (3) 



i=l 



Examples for basket weights are given in Table 1. For the DAX the weights are listed in 
Table 2. 

Theoretically the decomposition could change during the life of the option. However these 
events are relatively rare so that we ignore potential effects for the pricing in this paper. Let 
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Ticker Wgt 


Ticker Wgt 


Ticker Wgt 


ADSG.DE 2.20 
BMWG.DE 3.66 
CBKG.DE 8.20 
DBlGn.DE 2.12 
EONGn.DE 18.76 
HNRGn.DE 0.68 
LING.DE 1.92 
MRCG.DE 0.74 
RWEG.DE 4.71 
SIEGn.DE 9.28 


ALVG.DE 5.16 
BAYG.DE 9.41 
DAIGn.DE 9.72 

DPWGn.DE 9.57 
FMEG.DE 2.14 

HNKGp.DE 2.03 
LHAG.DE 5.21 
MEOG.DE 1.26 
SZGG.DE 0.43 
TKAG.DE 3.79 


BASF.DE 10.46 
BEIG.DE 0.93 
DBKGn.DE 7.04 
DTEGn.DE 33.91 
FREGp.DE 0.92 
SDFG.DE 1.40 
MANG.DE 1.13 
MUVGn.DE 2.25 
SAPG.DE 9.63 
VOWG.DE 0.98 



Table 2: Decomposition of DAX on the 31 July 2009. The first column displays the (Reuters) 
names of the assets followed by the weight in the index. 



us denote by Ci{K,T) the options prices of the i-th constituent. Similar to before, these 
prices are observed in the market so that individual distribution can be extracted according 
to 

^W)^.*,...-). (4) 
The integration of Eqj2] and Eq|4] yields for the cumulative distributions 

Let u denote a n-dimensional standard normal correlated Gaussian variate with 

E^^iu,) = (7) 
E^iuf) = 1 

for a given correlation matrix pij. If one assumes a Gaussian copula one can construct the 
joint distribution from the constituents by inversion of 

n, = Niu,) = PiiSiiT) <Ki) t = l,...,n. (8) 

The resulting multivariate distribution is denoted by 

pro6(T,5i(T),52(T),...,5„(T)). 

Hence, index option prices can be computed in two different ways: 

1 . by means of Eqjl] representing the market quotes of the index option 
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Figure 1: Comparison of implied volatilities for the two year option market of EuroStoxx 50 

between the market quotes and option prices obtained via multidimensional integration 
assuming a Gaussian copula. Only (roughly) half the skew is accounted for by the skew 
of the individual distributions. 



2. by means of multidimensional integration over its constituents, e.g. 

_ poo f " 1 ^ 

C{K, T) = DF y dSi...dSn prob{T, S,, S2, . . . , Sn) |5] "^'^'^ " ^| (9) 

If a Gaussian copula was to represent an adequate picture for combining individual spot 
distributions one should find 

C{K,T) = C{K,T) yK,T. (10) 

We have carried out this calculation, with some of the results depicted in Fig.l, Fig.2, 
Fig. 3. However instead of comparing option prices or distributions directly it is more mean- 
ingful to compare implied volatilities instead. 

Fig. 1 shows that individual option distributions fail to explain the volatility skew of the 
index for both EuroStoxx 50. Only roughly half of the skew can be attributed to the skew 
of the individual assets themselves. Note that this inconsistency does not necessarily imply 
dispersion arbitrage opportunities between the index option market and the options markets 
of and the individual assets. It is more reasonable to conclude that the option market data 
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Four-year Implied Volatility for .STOXX50E 
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Figure 2: Comparison of implied volatilities for the 4-year option market of EuroStoxx 50 and 
results obtained by integration of individual distributions that were combined by a 
Gaussian copula. 



suggests a type of dependency between individual stocks that is more complex than the one 
described by a Gaussian copula. 

In order to study whether the effect is persistent in time Fig. 2 displays the result of 
the 4-year market. The qualitative picture remains the same: Deviations from normality of 
individual log-returns cannot explain the non-normality of log-returns of the index. 

The results presented so far could represent an anomaly of EuroStoxx 50. However Fig. 3 
displays the result of calculations for the DAX that are qualitatively similar. 

Fig. 1,2,3 show that individual option distributions fail to explain the volatility skew of the 
index for both EuroStoxx 50 and DAX. Only roughly half of the skew can be attributed to 
the skew of the individual assets themselves. Also, other choices for the correlation matrix 
in Eq|7] fail to generate enough skew. Possible improvements by "scanning" the space of 
correlation matrices would be marginal at best. Note that reasonable choices for correlation 
matrices should result in basket implied volatilities that match the index's implied volatility 
curve around the most liquid point, which typically occurs at a strike of about 100%. 

Thus, the option market data strongly suggests the behaviour of dependencies between 
individual stocks in the index are more complex than the one described by a Gaussian copula. 
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Two-year Implied Volatility for .GDAXI 




• Market Implied Volatility 

* Gaussian Copula Implied Volatility 
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Figure 3: Comparison of implied volatilities for the 4-year option market of DAX and result 
obtained by integration of individual distributions that were combined by a Gaussian 
copula. 
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3 One parameter family of correlation matrices 



In this chapter we lay down the groundwork for the local correlation model (LCM) described 
in the next chapter. We define a family of perturbations of a given (center) correlation matrix 
that preserves all the relevant properties of the system, such as its positive semi-definiteness. 
The construction of this family is motivated by Poisson jump processes (common jumps gen- 
erally increases the correlation between assets while un correlated jumps cause de-correlation) 
even though jumps do not play an explicit role in the LCM model. The important result of 
this chapter lies in the fact that a one-parameter family of correlations can be defined that 
allows one to continuously dial correlations in between the identity matrix and a correlation 
matrix that is identical to one, i.e.pij = 1. 

Definition 3.1. Let r^") be a n- dimensional vector space in R with a canonical basis given 
by {e-i : i = 1, . . . ,n} . A linear operator that is represented by a n x n matrix C in the basis 
{cj} is called a correlation matrix if 

(i) C is positive semi- definite, i.e. u^Cu > Vz/ G V^*-"-* 
(a) All diagonal elements of C are one , i. e. Ca — 1 i — 1, . . . , /i 

Remark 3.1. Since any projection of C onto the subspace ofV^""^ is positive semi-definite, 
for any pair i,j G {1, . . . , n} the 2x2 matrix M defined by 



is a correlation matrix as well. As detM > it follows that — 1 < Cjj = Mj^q^^j^q-^ < 1 
Hence from Definition 3J_ it follows that all off-diagonal elements of C must be between -1 
and 1. 

Lemma 3.1. Let be a n-dimensional vector space in R and p a correlation matrix . 
Let u G V^'^'' be an arbitrary vector and k G {0, 1}. The matrix defined by 

h = { (11) 



is a correlation matrix. 




Proof: In order to prove Lemma 3J^, it suffices to construct, for a fixed time T = 1, n 
random variables i = 1, n with 

< E'3[(X^ - E«(X^))2] < oo and E«[(X^ - E«(X^))(X^ - ^^(X^))] < oo 

such that 

E<^[{X^^-E<^{XmX'r-E<^{Xi.))] ^^^^ 
EQ[{X^ - E«(X^))2]EQ[(X^ - EQ{X^)y] 
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yields EqjTT] Let {^l,Ft,Q) be a probability space where the filtration Ft is generated by 
a n-dimensional correlated Brownian motion ut together with two Poisson process Ji = 
iVf Vi, J2 = A^f V2 of (the same) intensity A where A^^*^ i = 1, 2 is the number of jumps of 
Poisson process i that occurred between and t and z/j is the jump size. Ut is given by 



d < ujI.ujI >= Pij(t) dt 



(13) 



If one defines 



It is straightforward to show that 



E[Xt. 



-XTui 



and 



E[(X 



T + {vif\T{l + \T). 



In the case where both jumps are triggered simultaneously, e.g. A^^ 



(1) 



(2) 



computation of Eq. 12 yields 



In the case where both jumps are independent one finds 

pij 



a direct 



(14) 



Pij 



(15) 



Equations 14 and 15 can be combined with the result 

Pij ~\~ KiUiUj 



Pij 



(16) 



where k = {0, 1} and Ui = \f\vi □ 



Remark 3.2. Equations I4 and 15 are consistent with intuition: In cases where sgn{i'i) = 
sgn{vj), correlated jumps increase term correlations between two random variables whereas 
uncorrected jumps tend to decrease correlations (in absolute terms). 



Remark 3.3. One can reduce Equation 16_ to a 1-parameter family of correlations with 
parameter u & R by introducing a (fixed) "principal mode" ^ G V^'-"-' to describe perturbations 
Ui = u of the correlation matrix, e.g. 



PijW 




^ 7^ J e {1, . . . ,n} 



(17) 



This motivates the following 
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Definition 3.2. The set 

(18) 

where pa = 1 is called a one- dimensional family of correlation matrices with center p{u) = 
{p)ij hi = I7 (^iT'd mode ^. 



Remark 3.4. In the special case of a flat mode, i.e. ^ = (1, 1, . . . , 1)-^ Equation 18 reduces 
to 

.^jG{l,...,n} 
PiA^) ^ < (19) 
I 1 ^=3 

This particular form will prove very useful for the construction of the LCM model because 
of the analytic tractibility that goes along with it. This will turn out to be crucial if one 
wants to construct a model that will be of practical use. 
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4 Local Correlation model (LCM) 



The goal of this chapter is to extend the classic multi-asset Dupire local volatility model 
by introducing the notion of local correlations into the dynamics between the assets. Local 
correlations imply that, during an infinitesimal time-step, assets evolve according to a cor- 
relation matrix that depends on the current state of the system and generally varies from 
one time-step to another. The main result of this chapter is the derivation of a particular 
functional form of the correlation matrix that achieves consistency with the index skew by 
construction. This is crucial if one wants to avoid lengthy fitting procedures of correlation 
parameters. 



Assumption 4.1. Let {Q,Ft,Q) be a probability space where the filtration Ft is generated 
by a n- dimensional correlated Brownian motion Wt = {uj}, e.g. Ft = F^ where 



d < ujI,uj{ >= pij{t) dt (20) 



Pij(t) is a correlation matrix with Pijit) > 0. There exist n assets Si{t) i = l,...,n that 
evolve according to the following local volatility dynamics: 



-f = /ijdt + a,(t, SI) dul z = l,...,n (21) 
where p\ denotes a deterministic drift. The zero index defines a basket 

n 

S',^Y.^,Sl (22) 

together with 



dS^ = J2 (^idSi (23) 



i=l 



for some given basket weights < a* < oo. Locally there exists a Brownian motion ui^ 



that matches the dynamics Eq.23 in a weak sense, e.g 



-^ = p'tdt + ao{t,S^)du;'t (24) 



12 



Furthermore for each asset i = {0, ...,n} there exists a market that trades plain vanilla 
European options at all maturities T < T* < oo and strikes K < K* < oo. The set of call 
options with price Cl{T, K) at time t of the i-th asset is denoted by 



{Ci{T,K) : T<T*; K < K*; t = 0,...,n} (25) 



where 



dCi{T,K) d^CiiT,K) 

exist and are continuous. For a single economy {i = 0, ...,n) consisting of a determin- 
stic bond B{t,T) = exp{— (r^ ds)) together with SI , i.e. Si = {B{t,T), SI}, there exists 
an equivalent martingale measure (EMM) Qi. 

Lemma 4.1. Under assumption \4.1\ each single option market Eq.25_ has no arbitrage if the 
following statements are correct: 

(^) 

f^l = rt- ql (27) 
where rt denotes the (deterministic) instantaneous interest rate and ql the dividend yield. 



in 



(^i{t,S,)- , a^ciiT,K) i-0,...,n {21 

2^ dK'^ 



Lemma 4.1 states the Dupire-equation [7]. As there is some confusion in the hterature about 



the drift terms we repeat the proof in Appendix 1. 



Remark 4.1. The lack of continuity condition of type Eq.26 generally leads to serious im- 
pairments in the model's ability to calibrate to the market. Mostly this occurs when market 
participants define their own internal volatility surface representation that is based on param- 
eterizations that include sharp cut-offs or by kludging together parameterizations of different 
types at some point of the curve in a non-differentiable way. 

The following lemma deals with cross-arbitrage (also called dispersion-arbitrage ) between 
the individual option markets and the option market of the basket itself 

Lemma 4.2. Under assumption \4 1\ there is no-dispersion arbitrage between individual op- 
tion markets and the option market of the basket if and only if 



S',li',=Y,o^.-^l^tSl (29) 



i=l 
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and for any t <T* and for any [S], S"") G i?"'"'" one has 

^ = J2 «^ ^*^^(^' Si) ^jit^ Si) P^t) (30) 

, 1=1 / i,j=l 



As the basket weights are assumed greater or equal to zero in Eq,22, the maximal covari- 
ance is achieved in the case when the correlation matrix is identical one, i.e. pij = 1 . The 
lower bound is given in the case when the correlation matrix equals the identity matrix due 
to the non-negative correlations assumption made in Assumption 4.1 This gives rise to the 
following 



Corollary 4.1. Under assumption ^.i a necessary condition for no-dispersion arbitrage 
between individual option markets and the option market of the basket is 



^ {a,Sla,{t, Sl)Y < i J2 (^^Sl j alit, S^) < a. a,Sl S^t, S^) cr,{t, S{) 
1=1 \j=i / *)i=i 

(31) 

Note that EqjSl] could be violated by the market since all variables are infered from it. If 
this happens it would indicate arbitrage opportunities provided the form of Eq|2T]is correct. 

Proof: Note that from Ito's lemma one obtains 
(5° - Ky = C?{T, K) + dCl{T, K) 

= C?{T, K) + '-^^{dSl - f^^S^du) 

= C'uiT, K) + '-^^{dSl - f^^S^Ju) 
where d < S^, 5° >= du {ao{u, S^)S^)'^ which follows from Eq, 
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The cancellation of the last two terms is due to the Black-Scholes partial differential equation 
but can also be viewed as a direct consequence of the martingale representation theorem. 
The last line describes the hedge-replication strategy of the option payout in terms of delta- 
hedging with the basket spot. In order to prove the absence of dispersion arbitrage one needs 
to re-write this equation as hedge-integral in terms of its components itself. If one inserts 
Eq,23 one achieves this goal as long as the drift terms still cancel so that discounted option 
prices remain martingales. This is indeed the case as long as the quadratic variation of the 
basket obeys 

n 

dt {ao{t, S^)S^f = d< Sf, >= dt J2 «i (^jSi SicTiit, Si) aj{t, S{) pij{t) 
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Hence the above hedging strategy can be re-written in terms of an arbitrage-free rephca- 
tion strategy that hedges into the individual components instead: 



□ 



Theorem 4.1. (LCM) Let F^lpo,^) be a one- dimensional family of correlation matrices 
with center po and mode C, = (1, 1)^ and Pij{u) G F(po,0- -^^^ 



n 

covo = ^ aia,jSlSl(Ti{t, Sl)aj{t, Sl){po)ij (32) 



covi = aiajSlSlai{t, Sl)aj{t, SI) 



(33) 



Under assumption \4.1\ and the assumption of Eq\31j the option market of Eq\2^has no arbi- 
trage if Eq'2^^^E^29^ hold and the correlation pij in Eg. 30 is set to 



where 



= < 



covo-affi,S<^){S'>r 
■co^i-<72(t,5»)(5»)2 



if covo-al{t,S^,){S',f <Q 



else 



(34) 



K is set to one when covq — (7Q{t, S^){S^) < and zero else. 

Proof: The no-arbitrage hypothesis of individual markets Ei follows directly from Lemma 
4.1 In order to proof the absence of dispersion- arbitrage, because of Lemma [4]2 , it suffices 
to construct a correlation pij that satisfies Eq|30] explicitly. If one inserts Eq, 19 into Eq, 30 
one observes the following cases: If covq < (To(t, 5'P)(S'f )^ the center correlations are not large 
enough to satisfy Eq,30 and one needs to set k = 1 in EqjlQ] as pij 



1 for u — s> oo. Note 



that the covariance matrix defined by 



15 



ipo)ij + KU'^ 

1 + ^2 



(35) 



is a strictly increasing function in u"^ for k, = 1 and a solution must exist in this case due 
to the bound in EqjsT} The result is given by the first line of Eg. 34 



Similarly if the center correlation is too large implying that covq > a^lt, Sf){S^Y one 
needs to set k = in this case as pij — > for u — > cxo and i ^ j- cov^u^) is a strictly 
decreasing function and hence invertible in this area with a solution sure to exist due to 
the lower bound of EqjsTj The result is given by the second line of Eq,34 We note here 
that in the case where non-flat modes are defined (see Eq, 17) similar arguments apply for 
the existence of solutions to Eq|30] provided > 0. The explicit solution may require the 
application of a simple numerical root-finding procedure however □. 

Remark 4.2. For simplicity we restricted ourselves to a situation where pij > since this 
is generally fulfilled in the equity market. However assume that, instead of the bounds given 



by Eq.31, more general bounds are obeyed by the market 



^ a, a,Sl Sia,{t, Sl) a,{t, S{)f4^-- < ^ a^Sl al{t, S^) < «^ Sia,{t, SI) cr,{t, S}) 

i,j = l \j=l / *ij = l 

(36) 



where pf"^"" is a correlation matrix that may contain negative correlations entries. In 



order to accomodate this case also Eq. 11 needs to be generalized according to 



PijW 



up 



1 



^ 7^ J e {1, • ■ ■ 



3 



(37) 



This can easily be achieved by correlating not only the Bronian motions by p"-"'^'^ but also 



the changes in the counting processes dNi by means of p"^ in the proof of Lemma 3A 
The corresponding equation for u* are quite similar to before and given by 
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Two-year Volatility Skew of .STOXX50E in the LCM 
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Figure 4: The calculation compares the market implied volatility quotes for 2-year options with 
the (standard) multi-asset local volatility model as well as LCM. Whereas LCM tracks 
the market distribution quite closely, multi-asset local vol deviates from the market 
substantially. In order to identify possible sources of numerical errors, we included 
EuroStoxx 50 in the same simulation as a "single" asset as well. 



5 Simulation Results of LCM 

In chapter [2] we presented the results of the Gaussian copula calculation and demonstrated 
that the volatility skew of the components of an index do not suffice to explain the traded 
skew of the index itself. For the results presented here, the choice of correlation matrix was 
made similar to other market-participants even though the qualitative results are indepen- 
dent of the choice of correlations. 

In this chapter we show the results of a similar calculation using the multi-asset local 
volatihty model described in Eqj2T} when the basket consistency equation [30] is not enforced. 
As one can see in Fig. |4] the results are qualitatively similar to the copula calculation of 
chapter [2] A constant correlation cannot explain the steepness of the volatility skew of the 
index. 

In addition Fig. |4] presents results using the LCM model. LCM closely tracks the distri- 
bution of the index. 

LCM can provide information for the average (off-diagonal) correlations for a given op- 
tion. The results are presented in Table [3j The results show that for a two-year put option 
struck at 70% of its spot the correlation is 15% percentage points higher than for a corre- 
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Strike in terms of spot 


average correlation in percent 


0.7 


58.3 


0.8 


56.1 


0.9 


53.8 


1.0 


50.6 


1.1 


46.9 


1.2 


43.5 



Table 3: Average instantaneous correlations of two-year options between the constituents of Eu- 
roStoxx 50 along the simulated paths in LCM. Correlations are significantly higher on 
the downside. 



Strike in terms of spot 


average correlation in percent 


0.7 


55.0 


0.8 


51.2 


0.9 


47.7 


1 


44.5 


1.1 


41.9 


1.2 


39.5 



Table 4: Average instantaneous correlations of four- year options between the constituents of DAX 
along the simulated paths in LCM. Correlations are significantly higher on the downside. 



spending call option struck at 120% of spot. To our knowledge this strike-dependence of 
correlations is substantially higher than what a trading desk typically accounts for in taking 
reserves for correlation risk (typically 5% points for liquid instruments) . Hence, one expects 
mispricings in cases where an exotic payout is convex. The corresponding results for 4-year 
options are given in table 4. 
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6 Implementation details of LCM 

In this chapter we briefly outhne how Monte Carlo code needs to be structured to implement 
LCM. 

Consider a simulation of nPaths, nAssets and nDates. Traditionally Monte Carlo code 
generates one path at a time for each asset in turn. However this will not work in this case 
as the choice of correlation matrix depends on the state defined by all the assets. Hence the 
structure should look as follows: 



begin 

call PrepareCho\esky{numberOf States, shift, centerC orr elation) 
for ipath := 1 to nPaths 

currentPath = initialSpots 
for idate := 1 to nDates 

randoms = GenerateRandomNumbers() 

StateVar = CalculateCorrelationState(cMrrentPat/i) 

CorrelatedRandoms = CorrelateRandoms(ran(ioms, C holesky Decomp[StateV ar]) 
for iasset := 1 to nAssets 

PathArray[ipath] [iasset] [idate] = GeneratePathAnay{CorrelatedRandoms) 

end 

currentPath = PathArray[ipath] 

end 

end 

call CalcPayout(Pait/iylrra?/) 
end 

In order to correlate random numbers one needs to perform a Cholesky decomposition. 
In the case where correlations are constant the latter needs to be computed just once, prior 
to the simulation of the paths. This is not true for LCM as every set of paths will generally 
require a different matrix to be decomposed. In order not to adversely impact performance 
it is useful to prepare a table of correlation matrices beforehand together with their Cholesky 
decompositions. The table implements a discrete version of F{pbase, ^ = (1, 1)^) for pijiu*) 



from Eq. 19 for about 50-100 values for u* , along the lines of: 



begin PrepareCholesky{numberOf States, shift, center Correlation) : 
/ / This function prepares a list of an odd number of correlation 
/ / matrices together with its Cholesky decomposition and stores the 
/ / results together with the state variable. The index k = k* = (numberOfStates - 1) / 2 



/ / corresponds to centerCorrelation. For k > k* one sets k = 1 in Eq, 34 
/ / whereas for k < k* it is set to zero in this case, 
if number O f States = even 

numberO f States := number O f States + 1 
end 

m := {numberOfStates — l)/2 
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correl[m] = center Correlation 
for := m + 1 to number States 

I := k — m 

Sum Over ij 

\f i = J 
correl[i][j] = 1 

else 

u = shift + / * shift 

correl[k][i][j] = correl[k][j][i] = {centerCorrelation[i][j] +u^)/{l + u^) 
StateVariable[k] = u 

cholesky[k] = CholeskyDecompose{correl[k]) 

correl[m — l][i][j] = correl[m — l][j][i] = centerCorrelation[i][j]/ {1 + u"^) 
StateV ariable[m — I] = —u 

cholesky[m — I] = CholeskyDecompose{correl[m — I]) 
end 

end 
end 

Once the simulation has started the function CalculateCorrelationState{currentPath) 



simply calculates the state variable according to Eq, 34 and looks up the closest position in the 
pre-calculated table from PrepareC holesky . Note that all current path values together with 
each constituent's local volatility and the local volatility of the index need to be available to 



compute Eq,34 
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7 Residual (exotic) Correlation Risks in the Pricing of 
Multi-asset Derivatives: The "Chewing- Gum" Ef- 
fect 

LCM allows one to compute single stock delta and vega exposures for multi-asset derivatives. 
In addition it produces sensitivities when the basket volatility skew and thus the distribution 
of the basket is changing. As the basket volatility depends on the correlation, vega hedging 
of the index component of an exotic, to some extent, hedges the correlation exposure as well. 

In a sense the idea is to "lock in" the projection of the exotic instrument onto the bas- 
ket "subspace" and therefore to hedge some of the correlation exposure along with it. 

However a hedge into the index only hedges overall (correlation) movements of the market. 
Mathematically, the specification of the distributions of the sum together with its components 
does not suffice to fully describe the joint distribution. We give a simple example in footnote 

i- 

In order to demonstrate this point, consider a basket Bt of two equally weighted stocks 
Si{t),S2{t) with Si{0) = 100,52(0) = 100 , e.g. 

Bt = liS^it) + S2{t)) 

and consider a "chewing-gum" move where Si{t) goes to 80 and 5*2 (t) to 120. As the 
basket level stays the same in such moves, neither the individual distributions nor the basket 
distribution change in this case. However an option that pays the worst of the two assets is 
affected in this case. 



In this chapter we try to get an estimate of this residual risk by finding alternative 
solutions that provide consistency with the basket skew. Any price sensitivity of worst-of 
options will be attributed to exotic correlation risk in this case. 

The LCM model was introduced for a family of correlations pij G F{p, ^) centered around 
p with a fiat mode ^ = (1, 1)^. The merits of this choice lie in the fact that p{u) can easily 



be inverted analytically in this case with the result given by Eq, 34 However the model can 



easily be extended to non-fiat perturbation provided > = 1, ...,n . This is feasible as 
long as one is willing to deal with a slightly more involved condition for the inversion of p{u) 
(see remark at the end of the proof of theorem 4.1). 

In this paper we pursue a simpler approach by simply varying the center-correlation instead. 
Recall that the center-correlation is merely a starting point around which a table of corre- 
lation matrices is constructed together with their Cholesky decompositions. The individual 



selection of correlation matrices is done locally according to Eq. 34 



In the following we present the results for two-year put options on the worst-of all 50 assets 
of EuroStoxx 50, which pays at maturity 

'^''''^5i(0)' -'550(0)^^ 
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Strike in terms of spot 


worst-of put price: setl 


worst-of put price: set2 


0.6 


0.306 


0.463 


0.7 


0.391 


0.560 


0.8 


0.482 


0.657 


0.9 


0.575 


0.754 



Table 5: Price of worst-of put 



for two different choices for the center-correlation: 

Set 1 presents results based on (1, 1)-^) where p was chosen close to where we 
believe the "market" was trading whereas set 2 is based on F((iiagf(l,...,l),(l,...,l)^) with 
the identity matrix chosen as center. 

We have also computed the average correlation which comes out roughly the same in 
both cases. Henceforth despite the fact that all individual as well as basket distributions 
are the same, the value of worst-of puts substantially differ between the two sets. In a 
sense basket options are mainly sensitive to overall shifts of the correlation matrix only and 
hence to the first principal component of the correlation matrix itself. This is different for 
worst (best)- of call and put options which show significant sensitivity to "higher" vibrations 
such as chewing-gum moves representing higher principal components of the correlation 
matrix. This is why liquid prices for worst (best) -of options would be very useful in providing 
further important clues about the joint dynamics of the system. In other words, the volatility 
skew for the basket can provide only little insight into the pricing of the worst-of (best-of) 
options, particularly when the number of basket members is large. This makes them ideal 
candidates to augment the correlation information provided by the basket. 
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8 Conclusion 



In this paper we show that the index option's market is inconsistent with the option market 
of the individual constituents if the dependency between the assets was to be described by 
a deterministic correlation matrix only or by a Gaussian copula. The data strongly suggests 
higher correlations on the downside and gives rise to a correlation skew which is computed 
explicitly in this paper. Roughly only half of the skew of the index can be attributed to the 
skew of the individual components, indicating that state-dependent correlation dynamics 
could play a crucial role in explaining the distribution of the index. 

This paper generalizes Dupire's local volatility model to the multi-asset case in a numeri- 
cally efficient way and provides consistency not only with the individual option's markets but 
also with the volatility market of the index itself by making correlations a dynamic variable. 

We also discuss the non-uniqueness of the solution and introduce the '"chewing-gum"' 
effect in this context where single stocks move against each other while leaving the distri- 
bution of the basket as well as all individual distributions unchanged. However worst-of 
options show price sensitivity in such moves which makes them ideal candidates to further 
"complete" the market and augment our knowledge of the joint dynamics of the system. 
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9 Appendix 1 



In the first step consider n + 1 separate economies Ei for i = 0, eacli one consisting of 
a deterministic bond B(t, T) = exp{— (r^ ds)) together with SI , e.g. Ei = {B(t, T), SI}. 
Because of the specific form of Eq|2T| that is the fact that volatihty depends on the spot 
only, each Ei is complete. The existence of an equivalent martingale measure states that the 
dividend- adjusted processes SI = Sl/e~'^^ are martingales under the cash numeraire B[t,T) 



yielding Eq,27 [6J. The specific form of Eq,28 has been derived in the literature already (see 



for example [7J). However this paper ignores drift terms. In several instances Eq,28 has 



been published incorrectly. The second term in the numerator ql Cl{T,K) is mistakenly 
replaced by /ij K). We quickly re-derive this equation at the end of this Appendix. In 

the second step one considers the economy e = {B(t,T), S^, , S^} that trades the option 
market specified in Eq,25 In order to avoid "spot" arbitrage between the index and its 



components taking appropriate {Qi ) expectation values on both sides of Eq j22| immediately 
yields Eq 29 One is left to show that Eq,30 avoids dispersion arbitrage between the basket 
option market and the individual option markets. This is done in Chapter |4] . 



Re-derivation of Dupire's formula: 
If one defines 



6 = {St - Kypt 

where Pt = exp{— j^Tgds). The application of Tanaka's formula yields 



dit = A [ls,>KdSt + -5{St - KySfaidt - n{St - K)' 
Stochastic integration on both sides with subsequent expectation taking yields: 



C^{T,K)-{St-K)- 

Differentiation on both sides with respect to T gives: 
dC,{T,K) 



dt (3t ^ltE[Stls,>K] + ^K^E[5{St - K)al] - nE[St - K)- 



dT 



(3t (/Ut E[STlsr>K] + \k^E[6{St - ir)4] - ttEISt - K^]^ 



Note that plain vanilla options obey 

Co{T,K) = (3T{E[STlSr>K] - KE[ls,yK]) 

The second expectation is just a call spread, e.g. E[1sj,>k] = — ^^°q^^^ ■ Similarly the 
expectation of the Dirac delta function can be inferred from option prices according to 
E[5[St — K)^] = ^ '^°2^^'^^ - If one inserts this back into the above equation and uses the fact 
that 

E[6{St - = E[6{St - K)]E\al\\ST = K] 



this yields Eq,28, as for a state-dependent local vol the expectation drops, e.g. E[a^\\ST 
K] = aUSr = K) 
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